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SUMMARY 


A  simplified  procedure  has  been  developed  to  calculate  the  1-P 
oscillating  aerodynamic  thrust  loads  on  single -rotation  propellers  in 
pitch  at  zero  yaw  on  tractor  airplanes.  The  application  of  this  procedure 
requires  only  a  knowledge  of  the  upflow  angles  at  the  horizontal  center 
line  of  the  propeller  disk.  These  angles  may  be  obtained  by  computational 
methods  thereby  obviating  the  need  for  experimental  flow  surveys. 

The  1-P  thrust  loads  and  1-P  vibratory  stresses  computed  by  the 
simplified  procedure  are  compared  with  those  calculated  using  all  the 
flow-field  characteristics  measured  for  several  wing-nacelle-fuselage 
combinations.  These  comparisons  show  good  agreement  for  the  1-P  thrust 
loads  and  the  1-P  vibratory  stresses. 

An  analysis  is  included  to  show  the  effects  of  the  various  flow-field 
parameters  on  the  1-P  thrust  loads. 


INTRODUCTION 


An  important  consideration  in  the  design  of  propellers  is  that  which 
deals  with  vibratory  stresses  which  occur  one  cycle  per  revolution  (com¬ 
monly  referred  to  as  1-P).  These  stresses  are  due  to  1-P  oscillating 
aerodynamic  thrust  loads  imposed  on  propeller  blades  as  a  result  of  pro¬ 
peller  thrust-axis  inclination  and/or  asymmetries  of  the  flow  fields  in 
which  the  propellers  operate.  Currently,  the  structural  design  criterion 
for  propeller  blades  is  based  upon  the  allowable  amplitudes  of  the  1-P 
vibratory  stresses.  Therefore,  there  is  need  for  a  simple,  rapid  method 
of  estimating  the  amplitude  of  the  1-P  thrust  loads. 

It  has  been  shown  in  previous  work  (refs.  1  and  2)  that  steady-state 
propeller  theory  may  be  applied  to  predict  adequately  the  oscillating 
portions  of  the  thrust  load  on  a  propeller-blade  element  around  the  disk 
if  all  the  characteristics  of  the  flow  field  in  which  it  is  operating 
are  known.  The  amplitude  of  the  1-P  component  of  this  variation  (the  1-P 
thrust  load)  then  may  be  obtained  and  this  load  used  to  determine  the  1-P 
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stresses  in  the  blade.  This  method  is  quite  laborious,  requiring  com¬ 
putations  for  at  least  1 6  points  around  the  disk  for  each  radial  position 
to  define  accurately  the  variation  of  the  oscillating  thrust  load. 

Further,  the  characteristics  of  the  entire  flow  field  at  present  may  be 
obtained  only  by  experimental  survey,  a  cumbersome  task  which  requires 
considerable  time. 

This  report  presents  a  procedure  for  calculating  the  1*P  thrust  load 
by  considering  only  the  upflow  angles  at  the  horizontal  center  line  of  the 
propeller  disk.  These  angles  may  be  computed  theoretically  by  the  methods 
presented  in  references  3>  and  5>  and  thus  the  need  for  experimental 
survey  is  obviated.  The  procedure  is  evaluated  by  comparisons  of  the  1*P 
thrust  loads  computed  by  the  procedure  with  those  computed  by  a  more 
rigorous  method  for  several  wing-nacelle-fuselage  combinations. 


NOTATION 


B  number  of  blades 

c  blade  chord 


blade-section  thrust  coefficient, 


section  thrust 
pn2D4 


ct 


l 


Actz 


D 

K 


blade-section  thrust  coefficient  due  to  section  lift  force  only, 
section  thrust  due  to  lift  (gee  Appendix  A. ) 
pn2D4 

incremental  values  of  the  coefficients,  the  value  at  any  angular 
position  minus  the  mean  value 

propeller  diameter 

Goldstein  correction  factor  for  a  finite  number  of  blades 


n  propeller  rotational  speed 

r  distance  along  any  radial  line  from  the  propeller  thrust  axis 

R  propeller  radius 

Vj  local  velocity  at  any  point  (r,fi)  at  the  propeller  plane 

(Direction  of  this  velocity  is  determined  by  the  angles 
9  and  ty.  See  fig.  1.) 


free-stream  velocity 

component  of  the  local  velocity  in  the  plane  perpendicular  to 
a  radial  line  (See  fig.  1.) 


velocity  ratio  at  any  point  (r,fl)  at  the  propeller  plane 

radial  location  of  any  blade  section,  ~ 

R 

the  fundamental  component  (one  cycle  per  propeller  revolution) 
of  the  oscillation  of  the  section  thrust  coefficient  or  of 
the  vibratory  blade  stresses 

geometric  angle  of  attack  of  the  airplane  thrust  axis 

propeller-induced  angle  of  inflow 

section  blade  angle 


tan 

amplitude  of  the  1-P  component  of  the  variation  of  incremental 
section  thrust  coefficient  due  to  section  lift  force,  further 
defined  as  the  1*P  thrust  load 

angle  of  outflow,  measured  from  a  line  parallel  to  the  thrust 
axis  in  a  plane  through  the  thrust  axis  (See  fig.  1.) 

angle  at  which  the  local  velocity  at  any  point  on  the  propeller 
disk  is  inclined  to  the  plane  perpendicular  to  the  radial 
line  through  that  point  (See  fig.  1.) 

mass  density  of  air  in  the  free  stream 

propeller -blade -element  solidity, 

27tr 

angle  of  rotational  flow  (an  apparent,  not  an  actual,  rotation) 
measured  from  a  line  parallel  to  the  thrust  axis  in  a  plane 
perpendicular  to  a  radial  line  (See  fig.  1.) 

(The  value  of  >| r  at  the  horizontal  center  line  of  the  pro¬ 
peller  disk  is  identical  to  the  upflow  angle  at  that  point 
and  corresponds  to  the  notation  "A”  in  ref.  6.) 

angular  position  about  the  thrust  axis,  measured  counterclockwise 
from  the  upper  vertical  position  as  seen  from  the  front 


-x 


blade-section  drag 
blade-section  lift 
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Subscripts 


t  effects  resulting  from  consideration  of  variations  in  i|r  only 

q}  V^/Vq  effects  resulting  from  consideration  of  variations  in  9  and 

V2/V o  only 

T  effects  resulting  from  consideration  of  variations  in  \|r ,  9, 

and  Vz/V0  acting  simultaneously 


PROCEDURE 

Development 


-Total 


Examinations  of  total  oscillating  thrust  loads  for  propellers 
operating  in  the  flow  fields  of  several  wing-fuselage-nacelle  combina¬ 
tions  have  indicated  h  simple  procedure  for  determining  the  amplitude  of 
the  1-P  thrust  load.  As  may  be  seen  in  the  following  sketch  of  a  typi¬ 
cal  total  oscillating  thrust  load  and  its  harmonic  components,  no 
significant  odd -order  sinusoidal  components  other  than  the  fundamental 
(l-P)  are  present  in  the  variation. i  The  lack  of  significant  higher 

odd-order  components  in  this 
variation  (and  in  all  other 
variations  examined)  makes  pos¬ 
sible  the  extraction  of  the  1*P 
^  component  of  the  thrust  load  by 

a  simple  operation  which  may  be 
applied  regardless  of  the  ampli¬ 
tudes  or  number  of  even-order 
components  present.  The  sim¬ 
plicity  of  this  operation  stems 
from  the  fact  that  all  even-order 
harmonics  repeat  their  magnitude 
and  sense  at  half-cycle  intervals 
of  the  fundamental;  whereas  odd- 
order  harmonics,  including  the 
fundamental,  repeat  their  magni¬ 
tude  but  are  of  opposite  sense  at 

half -cycle  intervals  of  the  fundamental.  Thus,  if  the  higher  odd-order 
harmonics  are  nonexistent,  one  half  the  arithmetical  difference  between 
the  values  of  the  total  oscillating  thrust  load  at  any  two  angular  posi¬ 
tions  l80°  apart  (i.e.,  half -cycle  intervals  of  the  fundamental)  will  yield 
the  exact  magnitude  of  the  1*P  component  at  those  points. 


90 

Angular 


180 

position,  D, 


270 

degrees 


^-Harmonic  analyses  of  some  of  the  variations  have  shown  third—  and 
fourth-order  components  of  small  amplitude  but  they  are  of  little  signif¬ 
icance  in  the  total  variation. _ _ _ _ _ 
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The  operation  may  be  simplified  by  utilizing  any  four  equally  spaced 
positions  around  the  propeller  disk.  For  the  specific  points,  0°,  90° , 
180  ,  and  270°,  the  expressions  shown  below  may  be  used  to  determine  the 
amplitude  and  phase  angle  of  the  1*P  component. 


1 

2 


£2=270°-! 


=  Sc-^cos  cpa 


(1) 


1 

2 


ft=o° 


M Q= 180°- 


=  Seisin  cpa 


(2) 


where 

the  amplitude  of  the  1-P  component  of  thrust  load 

Cpa  the  angle  between  the  position  of  the  maximum  1*P  magnitude  and 

the  90°  position. 

The  small  value  of  cpa  as  shown  in  the  foregoing  sketch  makes  plausible 
the  approximation 


&ctzc°s  cpa  «  &ctz 

hence,  only  the  value  of  c^.^  for  the  90°  and  270°  positions  would  be 

required.  The  procedure  then  consists  of  the  solution  of  equation  (l) 
rewritten  to  include  the  foregoing  approximation,  as  follows: 


ct' 


I2=9qO 


"  ct 


•l 


ft=270° J 


=  Set 


l 


(3) 


While  it  is  necessary  to  know  the  values  of  all  of  the  flow  parameters 
at  these  points  to  calculate  the  corresponding  values  of  ct  ,  certain 
characteristics  of  the  flow-field  parameters  make  possible  further  simpli¬ 
fications.  It  is  seen  in  the 
sketch  that  at  the  90°  and  270° 
positions,  Q  and  Vj/V, 0  are 
near  0  and  1.0,  respectively, 
while  is  near  its  minimum 
and  maximum  values.  Thus,  the 
values  of  Q  and  Vj/V 0  may  be 
approximated  as  0°  and  1.0, 
respectively.  Since  the  values 
of  \|r  at  these  positions 
may  be  computed  by  methods 
described  in  references  3  to  6, 
experimental  survey  of  the  flow 
field  would  be  unnecessary. 
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It  may  occur  to  the  reader  to  circumvent  the  approximation 

5ct.cos  cpa  *  bc±i 

by  utilizing  expression  (2)  and  the  flow-field  characteristics  at  the  0 
and  180°  positions.  However,  it  has  been  found  that  higher  odd-order 
components  which  are  insignificant  in  the  total  c-^  variation  did  not 

greatly  affect  the  accuracy  of  expression  (l),  but  their  presence  did 
greatly  affect  the  accuracy  of  expression  (2)  for  the  small  phase  angles 
cp  involved.  The  simultaneous  solution  of  expressions  (l)  and  (2)  with 
complete  flow-field  data  did  not  give  greater  accuracy  than  the  procedure 
to  be  examined  in  this  report. 


Evaluation 


As  a  means  of  indicating  the  accuracy  and  general  applicability  of 
the  procedure,  1-P  thrust  loads  computed  by  the  use  of  this  procedure  were 
compared  with  those  determined  by  a  l6-point  solution.  Measured  flow-field 
data  were  used  for  the  l6-point  solution.2  For  the  procedure,  measured 
values  of  \|r  were  used  and  the  values  of  0  and  V^/VQ  were  assumed  to 
be  0°  and  1.0,  respectively.  The  blade  stresses  (due  to  aerodynamic  load¬ 
ing  alone)  resulting  from  the  1-P  thrust  loads  computed  by  the  two  methods 
were  also  compared.  The  geometries  of  the  model  configurations  considered 
in  these  comparisons  are  shown  in  figure  2.  For  purposes  of  this  evalua- 
tion,  it  was  assumed  that  a  four-bladed  right-hand  propeller  was  installed 
on  each  model.  The  characteristics  of  this  propeller  and  the  method 
employed  in  the  computation  of  the  stresses  are  given  in  reference  1. 

iphe  i.p  thrust  loads  for  the  0. 7~i’adial  station  of  the  propeller, 
blade  for  all  models  as  computed  by  the  two  methods  are  compared  in  fig¬ 
ure  3.  Comparisons  of  the  radial  variations  of  the  1-P  thrust  loads 
computed  by  both  methods  for  two  models  are  shown  in  figure  4.  The  1-P 
blade  stresses  for  these  two  models  are  shown  in  figure  5.  From  these 
figures,  it  is  seen  that  the  procedure  is  generally  satisfactory. 


DISCUSSION 


Although  the  procedure  of  this  report  may  be  considered  sufficiently 
accurate  for  the  configurations  investigated,  it  was  concluded  that  i 


values  of  the  flow  parameters  used  in  these  computations  were 
determined  by  experimental  survey  and  were  reported  in  references  2,  3, 
and  6.  Since  data  were  obtained  at  the  port  propeller  planes,  the 
fields  for  the  starboard  propellers  were  assumed  to  be  the  reflected 
images  of  those  for  the  port  propeller. _ _ _ 
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would  be  impossible  from  a  cursory  examination  of  a  configuration  to 
predict  whether  this  accuracy  could  always  be  realized.  To  gain  an 
insight  into  the  cause  or  causes  of  the  1*F  thrust  loads,  it  was  decided 
to  investigate  the  separate  effects  of  the  various  flow-field  parameters. 
Therefore,  the  variations  of  the  total  thrust  loads  on  a  blade  element 
were  calculated  for  the  models  in  the  following  manner:  (l)  using  meas¬ 
ured  values  of  \|r,  and  assuming  0  =  0°  and  V^/Vq  =  1.0,  (2)  using'  meas¬ 
ured  values  of  0  and  V^/V o  and  assuming  \|r  =  0°,  (3)  using  measured 
values  of  \|/,  0,  and  Vi/Vo .  From  these  variations,  the  1-P  components 
were  extracted.  As  may  be  seen  from  figure  1,  the  first  case  shows  the 
effect  of  the  angular  variation  of  the  inflow  velocity  V1,  while  the 
second  case  shows  the  effect  of  the  magnitude  variation  of  the  inflow 
velocity  Vf .  While  it  is  not  apparent  from  the  equations  for  the  total 
thrust  variations  (see  Appendix  A)  that  this  method  would  yield  the 
exact  effects  of  the  variables,  it  was  found  that  the  summation  of  the 
thrust  loads  for  cases  (l)  and  (2)  was  nearly  identical  with  the  thrust 
variation  calculated  for  case  (3).  The  1-P  components  of  thrust  varia¬ 
tions  and  the  associated  flow-field  parameters  are  shown  in  figure  6. 

From  these  figures,  it  is  seen  that  \|r  produced  a  1*P  component 
which  peaked  at  approximately  90°  and  270°.  Since  the  procedure  of  this 
report  utilizes  only  \| r  at  the  90°  and  270°  positions,  it  essentially 
determines  the  1-P  thrust  load  due  to  From  these  figures,  it  is  also 

seen  that  0  and  V^/VQ  produce  a  significant  1-P  component  which  peaked  at 
approximately  0°  and  l80°.  The  models  which  have  a  large  1*P  component  due 
to  0  and  V^/Vo  are  characterized  by  a  relatively  large  blunt  nacelle. 
Further,  when  the  nacelle  is  mounted  on  the  wing,  the  magnitudes  of  these 
1-P  components  are  increased  due  to  the  velocity  gradient  induced  by  the 
wing.  It  should  be  noted,  however,  that  all  of  these  1-P  components  due 
to  0  and  V^/Vq  are  so  phased  as  to  add  little  to  the  total  1-P  com¬ 
ponent.  The  total  1-P  component  appears  to  be  more  dependent  on  the  phase 
angle  between  the  1*P  components  due  to  \|r  and  0  and  V^/Vq  than  on  the 
magnitude  of  the  1-P  component  due  to  0  and  V^/VG. 

Although  the  foregoing  analysis  has  shown  the  causes  of  the  1-P 
thrust  load,  the  inability  to  predict  the  phase  angles  involved  makes  it 
impractical  to  use  this  analysis  to  attempt  refinement  of  the  basic  pro¬ 
cedure  of  this  report.  The  analysis  may  prove  useful  for  estimation  of 
probable  accuracies  of  the  basic  procedure  for  other  configurations. 


CONCLUDING  REMARKS 


A  procedure  has  been  presented  for  calculating  the  1-P  thrust  load 
on  propellers  of  tractor  airplanes  in  pitch  at  zero  yaw,  whose  applica¬ 
tion  requires  only  a  knowledge  of  the  upflow  angles  at  the  horizontal 
center  line  of  the  propeller  disk.  The  evaluation  of  the  procedure 
showed  that  the  1-P  thrust  loads  and  the  1-P  blade  stresses  resulting 
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therefrom  could  he  predicted  for  the  airplane  configurations  investigated 
and  it  is  believed  satisfactory  results  could  be  obtained  for  similar 
configurations . 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Mar.  21,  1955 
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APPENDIX  A 


The  l'P  stress  problem  is  primarily  flatwise  bending  of  the  blade 
which  results  from  oscillations  in  the  blade  normal  force.  Therefore,  it 
is  the  sum  of  the  components  of  the  thrust  and  torque  forces  perpendicular 
to  the  chord  line  of  the  blade  element  which  must  be  determined.  The 
normal  forces  on  a  blade  element  may  be  expressed  as 

Nk  =  t  cos  p  +  fqsin  P  (Al) 

where 


t  section  thrust 

fq  section  torque  force 

The  section  thrust  on  a  blade  element  was  given  in  coefficient  form  in 
reference  2  as 


c  =  Kjt3x3  gi  (cot  q)  -  tan  y)  _  V8sin  t|t 


57-3 


,  0 W 

cot  *  *  — 


TtnDx 


(A2) 


where 


cp  =  cp_  +  ai 


cpQ  =  tan~xi 


V?cos  \|r 


itnDx  -  V?sin  \|r 


This  expression  may  be  simplified  to  consider  only  the  section  thrust 
due  to  lift. 


ctz  =  k*3*2 


57-3 


(cot  <p  +  ^ 


Vf  sin  \[r 
rtnDx 


(A3) 


This  is  the  expression  which  was  used  to  calculate  the  values  of  c^ 
considering  all  the  flow-field  parameters  as  measured  values.  The  values 
of  0°  and  1.0  may  he  substituted  for  6  and  V^/V o  in  this  expression  as 
described  in  the  text.  The  resulting  expression  is 
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where 


CtZ 


Kit3x3 


g-j 

57.3 


cot  cp 


cot  cp  + 


°1 

57-3 


V0sin  A2 
itnDx  J 


cp  =  cp0  +  ai 

,/  V0cos  \ 

cp  =  tan"1! — - - - — . - -) 

\jtnDx  -  V0sm  i| tj 


If  the  hlade  normal  force  is  expressed  in  coefficient  form  and  both  sec¬ 
tion  thrust  and  section  torque  force  are  expressed  in  terms  of  section 
thrust  due  to  lift,  cp,  and  7,  equation  (Al)  becomes 


c»b '  -  isr|)cos  p  +  cti(l  -  lr|)tan('p  +  r)sln  e  (A4) 

When  equation  (4)  is  rearranged 


_  CtZ  Jcos[ (7  -  cp)  -  2(3  -  cp)]  +  cos(7  +  cp) 


cos  3 


cos(7  -  cp)  +  cos(7  +  cp) 


(A5) 


For  simplicity,  the  term  in  the  brackets  was  assumed  to  be  unity  which 
is  tantamount  to  saying  that  the  drag  force  on  a  blade  element  does  not 
contribute  materially  to  the  blade-element  normal  force.  Therefore, 
only  the  thrust  due  to  lift  on  the  blade  element  was  considered  to  obtain 
the  normal  force  which  is  expressed  as 

CtZ 

CNb  *  cos  3 


V 
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( a )  Model  a 

Figure  2.-  Geometric  characteristics  of  the  models. 
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Figure  3.  —  Comparison  of  IP  loads  as  computed  by  the  procedure  of 
this  report  with  those  computed  using  the  entire  flow  field  for 
several  wing -fuselage-nacelle  configurations ;  x=0.7J  V0=l65mph, 
n=/250rpm,  /S^7R= 21.7 ° 
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Figure  5 —Comparison  of  IP  vibratory  stresses  resulting 
from  I  P  loads  computed  by  the  procedure  of  this  report 
with  those  resulting  from  IP  loads  computed  using  the 
entire  flow  field ;  Vo=!65mph,  n=/250rpm,  fl07R=2L7? acG=IO? 
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(b)  Model  c 
Figure  5—  Concluded. 
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Figure  6.  —  Concluded. 


